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ABSTRACT 

As pointed out in previous studies, the measurement of the skewness of the convergence 
field K will be useful in breaking the degeneracy among the cosmological parameters 
constrained from weak lensing observations. The combination of shot noise and finite 
survey volume implies that such a measurement is likely to be done in a range of in- 
termediate scales (0.5' to 20') where neither perturbation theory nor the hierarchical 
ansatz apply. Here we explore the behavior of the skewness of k at these intermedi- 
ate scales, based on results for the non-linear evolution of the mass bispectrum. We 
combined different ray-tracing simulations to test our predictions, and we find that 
our calculations describe accurately the transition from the weakly non-linear to the 
strongly non-linear regime. We show that the single lens-plane approximation remains 
accurate even in the non-linear regime, and we explicitly calculate the corrections 
to this approximation. We also discuss the prospects of measuring the skewness in 
upcoming weak lensing surveys. 

Key words: cosmology: gravitational lensing, large-scale structure of the universe 



1 INTRODUCTION 

Weak lensing by large scale structures has been recognized 
for a decade as a potentially powerful tool to probe the dis- 
tribution of dark matter in the Universe, as well as a way 



stressed that valuable constraints on cosmological parame- 
ters are expected to come from a joint measurement of vari- 
ance and high er order moments such as the skewness of the 
convergence ((Bernardeau, Van Waerbeke & Mellier 1997 
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Schneider et al. 1998 ; [Van Waerbeke, 
Bernardeau, Mellier, 1999])). This is because the variance is 



Bernardeau, Van Waerbeke fc McUicr 1997|; [Jain fc Seljak 



1997 ; gchncidcr et al. 1998[)). Reports of sound detections of 



s uch an effect have been done recently by differen t groups 
((|Van Waerbeke et al. 2000|; [Wittman et al . 2000|; |Bacon ' 



Refregier fc Ellis 200q ; |Kaiscr, Wilson fc Luppino 2000D ) 
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strongly dependent on both the amplitude of the mass power 
spectrum and the density of the Universe, while the skewness 
(properly normalized) is essentially a measure of the latter. 
Another motivation for measuring the skewness of the con- 
vergence comes from recent works on the possible existence 
of intrinsic correlations in the galaxy ellipticities ( ( 



Refregier fc Heymans 2000; Croft fc Metzler 2000; ICatelan 



ance of galaxy yllipLiciLies wiLli lespecL Lu a lauduni uiieuLa- 
tion of the intrinsic ellipticities, and are in very good agree- 
ment with theoretical expectations based on cluster normal- 
ized models of structure formation. These detections demon- 
strate that the systematic errors can be controlled down to 
a low enough level; this is a remarkable achievement that 
paves the way to genuine reconstruction of projected mass 
maps of the Universe. 

The aim of projected mass maps is not only to mea- 
sure the projected mass power spectrum. It has indeed been 
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lenge the cosmological interpretation of weak lensing surveys 
if part of the signal is due to intrinsic alignment of galax- 
ies. However, additional information on the statistics of the 
signal such as the skewness can provide a way to confirm 
the nature of it (e.g. a signal dominated by intrinsic ellip- 
ticity alignment s may not lead to a m easurable skewness as 
pointed out in ( |Croft fc Metzler 2000] )). 

Current skewness predictions rely mostly on perturba- 
tive analysis of structure growth which is expected to break 
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down at scales (around a few arcmins) where the lensing 
signal is easier to measure. It is therefore desirable to dis- 
pose also of intermediate and small scale predictions for the 
higher order moments. For the second moment, the variance, 
the nonlinear evolution of the powe r spectrum have been 
used since the earl y calculations (( Miralda-Escude 1991 
Jain & Seljak 1997; Schneider et al. 1998)). These calcu- 



lations are based on a fitting formula for the non-linear evo- 
lution of the mass pow er spectrum ((Hamilton et al. 1991 
Peacock fc Dodds 199q )) accurate to '^ 15% , dep ending on 
the cos molo gical model. (Jain & Seljak 1997) and (schneidoi 
et al. 1998) used this phenomenological approach to com- 



pute the non-linear convergence power spectrum, the small- 
scale shear variance and its correlation functions. The results 
are in remarkably good agreemen t with ray-tracing simula- 
tions ((Jain, Seljak & White 2000)). Unfortunately the non- 



linear evolution of higher order statistics cannot be simply 
described by the non-linear power spectrum alone because 
a description of the non- linear evolution of mode coupling is 
also required. For instance, the non-linear evolution of the 
skewness directly depends on how the bispectrum evolves. 

Prescriptions for the non-linear evolution of the bis- 
pectrum were not available until very recently. Such an is- 
sue is obviously closely related to the understanding of the 
non-linear gravitational dynamics that has been extensively 
studied for the last 10 years. Although these efforts were 
primarily done for studying the galaxy clustering properties 
one can take advantage of the results that have emerged 
in this domain. Recently, Hui (1999) extended the skewness 
calculations in the strongly non-linear regime using "hyper 



extend ed perturbation theory" (hereafter HEPT, (Scocci 



marro .'z Frieman 1999)). The predictions agree well with 
ray tracing simulations (( Jain, Seljak fc White 200(]| )) at 
very small scales (0.1'), but overestimate the simulation val- 
ues at larger scales. This is expected since HEPT was pro- 
posed to work in the non-linear regime, where the hierar- 
chical clustering ansatz is assumed to hold . Extensions of 
this approach presented in the literature (([Valageas 2000; 



Munshi & Coles 2000; Munshi fc Jain 2000a; |Munshi fc Jain 



2000b[) ^nd ( p3ernardeau fc Valageas 2000| )) are also based on 
the hierarchical ansatz, and thus suffer from the same limita- 
tions. Recently, Scoccimarro & Couchman (2000) obtained 
a fitting formula for the bispectrum based on numerical sim- 
ulations of CDM models. In this paper, we use this result 
to compute the skewness of the convergence k and compare 
it to ray-tracing simulations. This allows to bridge the gap 
between previous results restricted to either the non-linear 
regime or the weakly non-linear regime where perturbation 
theory holds. A different approach to this problem, using 
Press-Schechter halos, is given by Cooray fc Hu (2000). 

Taking advantage of these investigations, we also ad- 
dress the importance of the lens-lens coupling effects and of 
the Born approximation, which are routinely used in weak 
lensing calculations. It is usuall y assumed (see (Bernardeau, 
Van Wierbeke & Mellier 1997)) in semi-analytical predic- 
tions that the coupling between lens planes is negligible, as 
well as the successive deflections of the lensed light ray. This 
is equivalent to assume that the Universe can be approxi- 
mated by a single lens plane at some fiducial redshift. There 
is a priori no reason to believe that these approximations are 
still valid in the non-linear regime. Here we test quantita- 
tively these approximations down to very small scales both 



using our ray-tracing simulations and semianalytic predic- 
tions. 

The questions we address in this paper are therefore 
the following: does the non-linear prediction of the skewness 
S'3(k) based on the non-linear evolution of the bispectrum 
describe correctly intermediate scales? Is the skewness of 
the convergence still an efficient way of discriminating be- 
tween different cosmologies? Are the corrections to the un- 
perturbed ray-path given by the Born approximation and 
lens-lens coupling terms still negligible at non-linear scales? 

The paper is organized as follows. Basic definitions, per- 
turbative expressions and the description of the non-linear 
bispectrum are reviewed in Section 2. In Section 3, we give 
the relevant equations to calculate the variance and the 
skewness of the convergence, including correction terms due 
to lens-lens couplings and deviations from the Born approx- 
imation. Section 4 describes the ray-tracing simulations. We 
present our results in Section 5 and conclusions in Section 
6. 



2 DENSITY FIELD STATISTICS 
2.1 Perturbative Regime 

We assume the standard model of structure grow th from 



initia l Gaussian fluctuations. Perturbation theory ((Peeblei 



198C; Fry 1984)) which describes the evolution of the den- 
sity contrast (5(x, i) = {p{x,t) — p{t))/p{t) is well known, 
therefore we shall quote here only the results and notations 
needed in this paper. The Fourier transform 5{k,t) of the 
density contrast is defined as 



5(k,i) 



(27r 



(5(x, t) e ' ". 



(1) 



If (5 <c; 1 (perturbative regime), for an Einstein-de Sitter 
Universe (f2o = 1,A = 0), the solution of the equation of 
motion of 5 can be written as the following perturbation 
expansion 



5(k,i)=^a'(t)5«(k), 



(2) 



where a{t) is the scale factor such that a(to) = 1 today. 
For Q. ^ 1 the time dependen ce can still be facto rized to 
a very good accuracy (see eg. ( Bouchet et al. 1992 ) for full 
second-order results) although the growth factor is a more 
complicated function of time than a{t). At any order, the 
5'''(k)'s can be expressed in ter ms of the first order fluctu- 
ation <5(i'(k) (JGoroff et al. 1986|) 



5'''(k) = / d-'ki... / d^ki (5i5(k-ki-...-kO X 
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where the kernels _Fi(ki, ..., ki) are non-linear scalar func- 
tions of the wave vectors ki, ..., k^. The power spectrum P{k) 
and the bispectrum i3(ki,k2,k3) are defined as 

{5(ki)5(k2)) = fc(ki-fk2)P(fc) 

{5(ki)5(k2)5(k3)) = fc(ki-hk2 + k3) B(ki,k2,k3).(4) 



Incorporating Eq.(H) into Eq.(y) leads to the expression of 
the bispectrum as a function of the second-order kernel 
_F2(ki,k2) and the power spectrum 



reproduces HEPT for erg ~ 1. For more details see Scocci- 
marro & Couchman (2000). 
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Note that Eq.(W) implies that the bispectrum is defined only 
for closed triangles formed by the wave vectors ki,k2,k3. 

2.2 Non-Linear Power Spectrum and Bispectrum 

The non-linear prediction of the skewness requires two ingre- 
dients: the non-linear evolution of the power spectrum and 
bispectrum. To describe th e non-linear power spect rum, we 
use the fitting for mula of (([Peacock fc Dodds 1996 )), based 
on early work by ((Hamilton et al. 1991)), and obtained by 
fitting the non-linear power spectrum in numerical simula- 
tions. 

In order to describe the bispectrum at all scales, we 
use the fitting formula derived by Scoccimarro & Couch- 
man (2000) for the non-linear evolution of the bispectrum 
in numerical simulations of CDM models, e xtending previ- 
ous work for scale-free initial conditions (( gcoccimarro & 
Frieman 1999| )). The kernel F2(ki,k2) in Eq.(|) is simply 



3 GRAVITATIONAL LENSING STATISTICS 

Calculations of gravitational lensing statistics using pertur- 
bation theory have been presented in detail in the liter a- 
ture (e.g. see (Bernardeau, Van Waerbeke & Mellier 1997)). 
Therefore, we give without derivation the expressions of the 
second and third moments of the convergence field k as well 
as the skewness correction terms due to lens-lens coupling 
and deviations from the Born approximation. 



3.1 Basics of Gravitational Lensing 



We use notations similar to (Schneider et al. 1998) with 
c = Ho = 1, since the results do not depend on the Hubble 
constant. In a FRW Universe with a matter density Qo and 
cosmological constant A, the radial distance w{z) and the 
angular diameter distance fK{z) are defined as 
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replaced by an effective kernel F2 (ki, k2) such that 
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where Q.k = fio + ^a — 1 is the curvature. In the absence of 
lensing, the comoving angular distance x(0,w) at a radial 
distance w{z) between a light ray in the direction 6 and 
a fiducial light ray is x{6,w) = 6fK{w). In the presence of 
density fluctuations, a light ray experiences many deflections 
such that x(0, w) becomes 
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where $ is the 3-dimensional gravitational potential. The 
gradient V± denotes the derivatives with respect to x in 
the plane perpendicular to the line-of-sight at distance w. 
The amplification matrix A{6, w) describes the mapping be- 
tween the source plane at w{z) 



l + a^'-^iz)[0.7Qs{n)]'''{q/4) 

1 + (g/4)"+3-5 
1 + 0.4 (n + 3) g"+^ 

l_|_q,n + 3.5 

l + 4.5/[l.5 + {n + 3)^] (2g)"+3 



A{e 



9x 



(13) 



fK{w)de' 

which using the first order expansion of Eq.(h2h becomes 
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and q = k/k]siL{z), where 'i'Kk% i^Pi^{kisi l) ~ 1, and Phik) is 
the linear power spectrum at the desired redshift. The effec- 
tive spectral index is taken from the linear power spectrum 
as well. The function (53 (n) is given by 

(4-2") 
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The convergence k and the shear 'y are defined such that 
the amplification matrix takes the symmetric form 



which is the "saturation value" obtained in HEPT ((Scocci- 
marro fc Frieman 1999[)). From these expressions, it follows 



A{e,w) = 
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that at large scales, where the functions a = 6 = c=l, we 
recover the tree-level PT. On the other hand, at small scales, 
where a^ = (7/10)(53 (n)(T^ ' and 6 = c = the bispectrum 
becomes hierarchical with an amplitude that approximately 



K and "f are therefore line-of-sight integrals from the source 
distance w{z) to 0. If the sources are distributed in redshift 
according to some function ps{w), the convergence k can 
still be expressed as a simple line-of-sight integral 
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regime (see ( Bernardeau, Van Waerbeke fc Mellier 1997 ) and 



(pchneider et aL 1998)). The contribution to the skewness 



can be calculated by considering the perturbation expan- 
sion up to the second order in the gravitational potential in 
Eq.QlsL which leads to a second order correction term Al, 
in Eq.([l4|). After lengthy but straightforward calculations, 
the correction to the third moment {k'^)^™'^ can be derived: 



3.2 Variance and Skewness of the Smoothed 
Convergence 

In perturbation theory, the moments are calculated from a 
perturbative expansion of the density contrast. The i"^ order 
of the convergence field is therefore 
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The calculation of the smoothed second and third mo- 
ments is then straightforward. The leading order contribu- 
tions to these moments are {ii^)eo = {i^ )so and {k-'^)9o ~ 



3{k 



{1)^^(2) 



)gg respectively. Let us call I{k) the Fourier trans- 



form of the smoothing window and define the angular wave 
vector s = fK(w)li.±. As calculated in ( Bernardeau ^ Van 



second 'nioiiieiiL is 



Waerbepe fc Mellier 199V ) and ( gjhncider ct al. 1998D the 
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(20) 

These expressions rely on the validity of Eq.(|l4|) which in- 
volves two approximations: one is that we assume the light 
rays travel along the unperturbed light paths (the so-called 
Born approximation), and the second is that we neglect non- 
linear terms in the amplification matrix (which is physically 
equivalent to neglecting all the couplings between subse- 
quent lenses) . These two approximations ignore quadratic 
terms of the form k ■ k, k • 7 or 7 ■ 7 (see ( Bernardeau, Van 
Waerbeke & Mellier 1997 )) . There is no reason a priori to ne- 
glect these terms in the skewness calculation, because they 
give a contribution of the same order as Eq.(|20[). We shall 
address this question in the next Section. 

3.3 Skewness Corrections due to Lens-Lens 
Couplings and Deviations from Born 
Approximation 

The effect of dropping the lens-lens coupling and the Born 
approximation has been calculated in the weakly non-linear 
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The similarity with Eq.(M) is quite obvious, which shows 
that the two terms {{ii'^)0a and {k'')^™'^) are of the same 

order S^^' in the perturbative sense. The only difference 
between the two expressio ns lies in the double time inte- 
gral here. It was shown ((Bernardeau, Van Waerbeke & 



MeUier 1997| ) and ( pchneider et al. 1998] )) that this dou- 
ble integral makes the correction term Eq.(|2l[) only a few 
percent of Eq.(l2d). In other words, even though the two ex- 
pressions have indeed the same order ex S^^' , they differ by 
almost one order of magnitude just because of the additional 
fxiw — w')/fK{iu) factor in Eq.(|2l[) (in the case of a single 
source redshift this factor is equivalent to an additional g{w) 
factor, according to Eq.(O)). It is therefore correct to say 
that the lensing efficiency function g{w) can play the role 
of a perturbative parameter as the density contrast does 
(when S <C 1). The smallness of {K^)'e^" shows that even 
in the weakly non-linear regime, terms involving powers of 
g{w) higher than 1 will always be small compared to terms 
only proportional to glw). In the following we shall compute 
Eq.(bl|) in the non-linear regime, and use results from ray- 
tracing simulations in order to check if this statement is still 
valid when the density contrast is much larger than 1. 



3.4 Calculations and Filtering Schemes 

We can now calculate the variance, the skewness and its 
correction of the convergence field using the Eq.([19|), ( |20[ ) 
and (Ell) by forming the ratios: 
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It is worth to notice that AS'3(k) is weakly dependent on 
the cosmological parameters since the f2o dependence in the 
front of Eq.(|l9[) and Eq.fElh cancels out. The non-linear 
variance is given by Eq. (119) where the power spectrum is 
replaced by the Peacock & Dodds (1996) prescription as 
in Jain & Sel jak (1997). The non- linear third moment is 
given by Eq.(GG) where i*2(ki,k2) has been replaced by 
F| (ki,k2) and we use the non- linear power spectrum. For 
the skewness correction, we only need to use the non-linear 
power spectrum in Eq.(Ell). 



(21) 



We use two types of filtering schemes: top-hat filtering 
and Map statistics. The Fourier transform of the top-hat 
window is: 

Mv) 
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The Map statistic can be 



where J denotes a Bessel function 

defined either from the convergence k or from the tangential 

shear 7^: 



Ma 



d'e u{e)K{e) = / d?e Q{e)-it{e) 
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where e'^Q{e) = 2 J^ dO' 9' U{e')~e^U{e), and U^e) is the 
compensated filter. We use the compact compensated filter 
described in Schneider et al. (1998) where U{9) is defined 
only for 6 £ [0, 9o], and 
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ay 



The Fourier transform of U{6) is 



iiv) = 



12 Jijv) 



(25) 
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which is the function to be inserted in Eqs.(|l9|, EOL Ell). The 
motivation for using the compensated filters is illustrated 
on Fig. 2 of Schneider et al. (1998): this filter is a pass-band 
filter, which is a more direct estimate of the mass band 
powers than the top-hat filter. Another advantage of the 
compensated filter is that the measurement of the skewness 
does not require a mass reconstruction (whereas top-hat fil- 
tered statistic does) as it can be evaluated directly from the 
shear (Eq.(E4f)). Note also that the compensated filter de- 
fined above peaks at a scale ~ 0o/5, if Oo is the radius of the 
filter. Therefore at equal radius with the top-hat filter, Map 
statistics probes scales as small as ~ 6o/5. 



4 RAY-TRACING SIMULATIONS 

We use two different kind of ray-tracing simulations (high 
and low resolution simulations) which allow us to investigate 
scales from 0.1' to several tens of arcmins. The high resolu- 
tion simulations (0.1') were done by Jain, Seljak & White 
(2000), and are described in that paper. The investigated 
models in those simulations are 0CDM2 and rCDM (see 
Table bj). The area covered by each simulation is 9 square 
degrees. 

The numerical methods used for the low resolution sim- 
ulation are basically same as those of the high resolution 
simulation except for a use of a tiled set of the indepen- 
dent particle-mesh (PM) TV-body simulations^. The PM 
simulations use 256^ x 512 particles and a force mesh of 
the same size in a periodic rectangular comoving box, with 
outputs done along the light-cone (Hamana, Colombi and 
Suto 2000). We performed 10, 11 and 12 independent sim- 
ulations to generate the density field from 2 ~ 2 to the 



* This code, developed by the Numerical Investigations in Cos- 
mology group {N.I.C.), is an improved version of an earlier work 
by Bouchet, Adam & Pellat (1985) and by Moutarde et al. (1991), 
fully vectorized and optimized to run in parallel on several pro- 
cessors of a NEC-SX5. 



Table 1. Ray-tracing simulations used in this work. The SCDM, 
OCDMl and ACDM are described in Hamana et al. (in prepara- 
tion) and the OCDM2 and tCDM are from Jain, Seljak & White 
(2000). 9res indicates the angular resolution of the simulations 
Nrea indicates the number of realizations per model, and Area 
specifies the sky coverage for each realization in square degrees. 
Note that the Nrea = 40 realizations are not rigorously indepen- 
dent, being random rotations of a single suite of simulations. 



Simulations 


SCDM 


OCDMl 


ACDM 


OCDM2 


rCDM 


no 


1 


0.3 


0.3 


0.3 


1 


A 








0.7 








r 


0.5 


0.21 


0.21 


0.21 


0.21 


CTS 


0.6 


0.85 


0.9 


0.85 


0.6 


^s 


1.0025 


1.03 


1.0034 


1 


1 


Ores 


2' 


2' 


2' 


0.1' 


0.1' 


Nrea 


40 


40 


40 


7 


5 


Area 


25 


25 


25 


9 


9 



present for SCDM, OCDMl and ACDM models, respec- 
tively. We use the tiling technique first proposed by White 
& Hu (2000), i.e., the box size of each simulation is chosen 
so as to match to the convergence of the light ray bundle. 
The smaUest box size is comoving 80 x 80 x 160/i~^Mpc^ 
for three models, while the largest box size are comoving 
240 X 240 X 480/l-='Mpc^ 320 x 320 x 640/i-^Mpc^ and 
400 X 400 X 800/^.~^Mpc^ for SCDM, OCDMl and ACDM 
models, respectively. These tiled sets of sinrulations yield a 
field-of-view of 5 x 5 square degrees. 

The multiple lens-plane algorithm is used for ray-tracing 
as in Jain et al. (2000) . The lens planes are located between 
z = and z ~ 2 at intervals of comoving 80/i~^Mpc. We ob- 
tained 40 realizations by randomly shifting the simulation 
boxes. For each realization, 512^ rays are traced backward 
from the observer's point. The initial ray directions are set 
on 512^^ grids with a grid spacing of 5°/512 ~ 0.59 arcmin. 
Although the results of the ray-tracing simulations (the po- 
sition and amplification matrix of each rays) are stored at 
all lens planes, in this paper we only focus on the results at 
the plane closest to 2 = 1. The error bars on the results of 
the low resolution simulations correspond to the standard 
deviation among 40 realizations. 

The range of filtering scales where the results of the low 
resolution simulation are reliable is between 4' and 30' for 
the statistics of the top- hat filtered convergence with source 
redshift 2 ~ 1. On the other hand, for the Map statistics, 
that is between 20' and 100': the reason is that our com- 
pensated filter probes scales as small as ~ Oo/5 (see Section 
3.4), therefore Map predictions are likely to become inac- 
curate below 4' X 5 = 20'. The small-scale resolution limit 
comes from the finite spatial resolution of our PM simula- 
tions, while the large-scale limit comes from a lack of power 
in the density field on scales larger than the simulation box 
(finite size effects). 

In addition to the usual full ray-tracing simulations, 
we also performed approximated ray-tracing which neglects 
both the defiection of the ray trajectory and the couplings 
between subsequent lenses. To be more specific, the shear 
tensor, Um (see eq. 13 of Jain et al. 2000 for a definition) , is 
computed on the unperturbed ray position, and the distor- 
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Figure 1. r.m.s. (top panels) and skewness (bottom panels) of the convergence field as a function of top-hat smoothing scale for the 
OCDMl (far left panels) and SCDM (center-left panels) OCDM2 (center-right panels) and rCDM (far right panels) models given in 
Table ll. The solid line is the prediction using non-linear prescriptions for the power-spectrum and the bispectrum, and the dashed lines 
show the perturbation theory calculations. The vertical dotted lines in the left panels denote the reliable scale limit. The "large" error 
bars correspond to scaling the "small" error bars (thick lines) to a single 25 square degree (left panels) and 9 square degree (right panels) 
survey. 



tion tensor $„ in the right-hand-side of the eq. (14) of Jain 
et al. (2000) is replaced with the identity matrix to neglect 
all lens-lens couplings. The weak lensing statistics are then 
computed in the same manner as in the full ray-tracing case. 
The difference in the skewness obtained from the full and ap- 
proximated simulations are computed to test the validity of 
the perturbative treatment in the semi-analytic predictions 
(§5). 



5 RE JSULTS 



Figures |l| and Q show the measured r.m.s. and skewness of 
the convergence compared to the linear (dashed lines) and 
non-linear (solid lines) predictions. The square symbols rep- 



resent tlhe measurements done in the simulations. I'wo error 
bars are drawn per point: the small error bars correspond to 
the dispersion in the measurements taken into account all 
the realizations, whereas the large error bars are obtained 
by scaling (by \/Nrea) the small error bars to a single real- 
ization, relevant for assesing cosmic variance in a single field 
of the corresponding area. 

For the top-hat filtering scheme (Figure uh the non- 
linear predictions give very good results for all the scales 
investigated (i.e. down to scales as small as 0.1'). For the 



low rest)lution simulation (figure I lett panels), agreement 
starts tb ddtdtiot^te onCe wfi t^hCH the teliabfe fiCale lithit 
which is 4' (shown by the vertical dotted lines). At smaller 
scales the high resolution simulations results show that the 
fitting formula still works down to 0.1 — 0.2'. The agreement 



starts to deteriorate quickly above the scale of ~ 30' be- 
cause finite field effects become important, therefore these 
scales are not shown here. Moreover, from a practical point 
of view we are less interested in these scales because resid- 
ual systematics from shear measurements are supposed to 
dominate the lensing signal (Erben et al. 2000). The few 
percent discrepancy for the r.m.s. for the rCDM model (top 
right of Figure hi right panels) was already observed by Jain, 
Seljak & White (2000), but this is consistent with the ac- 
curacy of th e power spectrum fitting formula ((Peacock & 



Dodds 1996 )). The non- linear predictions for the skewness 
in Fig. Ill although consistent within the errors with ray- 
tracing simulations, shows fluctuations of S3{k) by approx- 
imately 10%. This is within the e xpected accuracy of the 



bispectrum fitting formula (15%, (Scoccimarro fc Couch 



man 200C)); however, one must also keep in mind that this 



formula was obtained from single realizations of CDM mod- 
els, and the bispectrum amplitude is rather sensitive to the 
presence of massive clusters (see Cooray & Hu 2000 for the 
effects of this on S'3(k)), therefore the residual 10% fluctua- 
tion might be due to cosmic variance. 

For the compensated filter (Figure bl), the difference 
between linear and non-linear variance is more significant, 
since when the spectral index becomes n < —2 the vari- 



ance in a compensated field starts to decrease ((Schneidei 



et al. 1998)), and this happens at much larger scales for the 



linear than the non-linear power spectrum. The agreement 
between predictions and simulations is very good for the 
low resolution simulations, including the skewness, and good 
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Figure 2. Same as Figure HI but for tiie compensated filter. 



for 0CDM2, the high resolution case. On the other hand, 
the agreement is not as good for the high resolution rCDM 
model. The likely explanation is that the small discrepancy 
observed in that case for the top hat filter (Fig. hh is ampli- 
fied by the compensated filter, which acts like a derivative 
operator on the convergence field: as a pass-band filter, the 
compensated filter is more sensitive to any small discrep- 
ancy (in fc-bands) between the non-linear power-spectrum 
fitting formula and the ray-tracing simulations, whereas the 
top-hat filter integrates over large range of fc's and tends to 
mask such discrepancies. We should note that the effects of 
resolution also cause the variance to decrease at small scales 
(top left panels in Fig. 0), therefore, at scales smaller than 
the resolution the agreement between the skewness predic- 
tion and the measurement in the low resolution simulation 
(bottom left panels in Fig. |2|) is coincidental. 

Figure M shows the skewness correction for top-hat 
smoothing in SCDM (top panel) and OCDMl (bottom 
panel) models. The dashed line and the solid line show re- 
spectively the predictions from perturbation theory and the 
non-linear fitting formula, whereas points with error bars de- 
note measurements in the ray-tracing simulations. Surpris- 
ingly, even in the non-linear regime, the correction remains 
small and, in fact, the relative skewness correction is smaller 
than in the weakly non-linear regime! As explained in Sec- 
tion 3.3, this means that coupling terms introduce a non- 
linearity which is always smaller than the non-linearity pro- 
duced by very dense (non-linear) structures. In other words, 
even in the non- linear regime, the dominant lensing con- 
tribution comes from the single lens-plane approximation, 
which physically means that the probability for a given line- 
of-sight to pass through two very dense structures at very 
different resdhifts is very low. Similar results hold for the 
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Figure 4. Skewness correction for top-hat smoothing AS3{k) as 
defined in Eq.Q for the model SCDM (top panel) and OCDMl 
(bottom panel). Squares with error bars correspond to measure- 
ments in the ray-tracing simulations, dashed lines to the predic- 
tion of perturbation theory, and solid lines to predictions from 
the non-linear fitting formula. 



Afap statistics (not shown here) where, in fact, the correc- 
tion terms to the skewness are even smaller. 

Figurerashows the predicted skewness S3{k) in the SI— A 
plane. Here we compare small scale (1') non-linear calcula- 
tions (left panel) to larger scale (10') perturbation theory 
calculations (right panel) for top-hat smoothing. Note that 
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Figure 3. Top-hat (left panels) and compensated (right panels) statistics for the ACDM model (table 



the choice of a scale for the perturbative estimate is almost 
not relevant since it does depend very weakly on scale (Fig- 
ure ol) . This Figure shows that the cosmological dependence 
is preserved in the non-linear regime, and that S^^k) re- 
mains essentially an estimator of f2 (the degeneracy along 
A is strong). This result has to be compared with the Cos- 
mic Microwave Background and the supernovae constraints, 
which are sensitive to the curvature and ~ A respectively. 
This clearly motivates further lensing studies as an inde- 
pendent means of measuring the cosmological parameters. 
In fact the lensing constraints are similar to constraints de- 
rived from the galaxy catalogues or the velocity fields, al- 
though the physics is much simpler here, and weak lensing 
does not suffer from the galaxy bias problem. 

The skewness predictions in the non-linear regime are 
twice their value using perturbation theory at large scales. 
Therefore the cosmic variance and noise analysis done in Van 
Waerbeke, Bernardeau & Mellier (1999) underestimated the 
capability of lensing surveys by a significant amount, essen- 
tially because they used very low resolution and approxi- 
mate dynamics in their simulationsPl. In this work the au- 
thors shown that a 25 square degrees survey would measure 
the skewness at 9cr, in standard observational conditions, if 
fl = 0.3 (3(T if f2 = 1). Here, the signal amplitude is twice 
larger, therefore we would need ~ 3 — 4 times less sky cover- 
age to reach the same signal-to-noise ratio to measure Sz{k?). 

T In this work the authors used very simple lensing simulations 
where the mass distribution was modeled with second-order 2- 
dimensional lagrangian dynamics, which they had to smooth at 
small scale in order to avoid singularities. Because of this smooth- 
ing they obtained a skewness of the convergence close to its value 
in perturbation theory. 



However, the exact gain depends on the precise amount of 
cosmic variance and shot noise at this scale. Since the cos- 
mic variance is already given in Figure nl we can simply 
add the shot noise obtained in Van Waerbeke, Bernardeau 
& Mellier (1999) to the high resolution simulation results 
presented here. It turns out that the shot noise is roughly 
(JSj ~ 10 at 1 — 2 arcmin, therefore we can reach the signal- 
to- noise (7 — 9cr for OCDM) with only 9 square degrees with 
a single measurement at 1'. However we should keep in mind 
that our results for the error on the skewness in the ACDM 
model are consistent although slightly smaller than previous 
estimates from White & Hu (2000), possibly due to the fact 
that our realizations are not completely independent (being 
random samples of a single suite of simulations). 

Finally, Fig. o shows the cross-correlation coefficient rtj 
between the variance measured at angular scale 9i and 6j, 
for the ACDM model (similar results hold for the other cos- 
mological models and for the skewness measurements as 
well). These results were obtained by averaging over the 
40 realizations in our ensemble. The top panel presents re- 
sults for the top hat filter, whereas the bottom panel cor- 
responds to the compensated filter. We show three angular 
scales, 9i — 5.6,10,23.7 arcmin for the top hat case, and 
6i — 23.7, 42.2, 75 arcmin for the compensated filter. The 
top hat measurements of the variance show a large cross- 
correlation, in fact, the cross-correlation coefficient is almost 
unity among all scales. On the other hand, as expected from 
being a band-pass filter, the compensated filter measure- 
ments are much less correlated. Note however than this re- 
sult only includes cosmic variance contributions, i.e. it does 
not include shot noise, which can be a significant source of 
cross-correlations at small scales. 





Figure 5. The expected skewness S'3(k) in tlie non-linear regime at a scale of 1' (left) compared to 53 (k) computed using perturbation 
theory at a scale of 10' (right) in the fl-A plane. Top-Hat filtering was used. The sources are located at rcdshift unity. Note that for a 
given O-A the skewness in the left panel is roughly twice the value on the right panel. 
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Figure 6. The cross-correlation coefficient rij between the vari- 
ance measured at angular scale 6i and 6j , for the ACDM model. 
The top panel shows rij for top-hat filtering at 8i = 5.6, 10, 23.7 
arcmin as a function of 9j . The bottom panel shows r^j for com- 
pensated filters at 6, = 23.7, 42.2, 75 arcmin as a function of 6j. 



6 CONCLUSIONS 

We calculated semi-analytically the r.m.s. and the skewness 
of the convergence field k and showed that the non-linear 
prescription for the 3- dimensional bispectrum given in Scoc- 

rimarrr) kr riniipliman (9000) Hparrihpa thp trangitinn frnm 



linear scales to the non-linear regime is clearly visible and 
also accurately described by our semi-analytical approach 
for both top-hat and compensated smoothing schemes for 
the OCDM, SCDM and ACDM models (The rCDM pre- 
dictions are not very accurate for the compensated filter, 
but we suggest that it might come from innacuracies in the 
non- linear power spectrum fitting formula). 

Our results extend previous non-linear calculations of 
the skewness of k which eith er assumed the validity of HEPT 
at all redshifts (( Hui 199£)) or re l ied on the hierarchica l 
clustering ansatz ((Valageas 2000 ; Munshi & Coles 200C; 
Bernardeau & Valageas 2000|)) at intermediate scales. Our 



semi-analytical approach allows us to investigate the cosmo- 
logical dependence of the skewness of k in as many models as 
desired, without having to perform a ray tracing simulation 
each time. 

By comparing the skewness in ray-tracing simulations 
with and without lens-lens coupling and Born approxima- 
tion terms, we demonstrated that the single lens-plane ap- 
proximation remains valid even in the highly non-linear 
regime, which is a non-trivial statement about weak lens- 
ing at small scales. We found that the contribution of these 
terms to the skewness is in remarkable agreement with the 
semi-analytical predictions. 

We have also calculated 5*3 (k) in the Q, — A plane and 
shown that the signal is about twice the expected value from 
perturbation theory with roughly the same degeneracy in 
A as in the perturbative case, generalizing previous results 
(Jain et al. 2000, White & Hu 2000). In addition, our re- 
sults on the cosmic variance of the skewness suggests that 
a combined analysis of the existing weak lensing surveys 
(dVan Waerbeke et al. 2000|; Iwittman et al. 2000|; [Bacon, 



perturh lative to non-linear angular scales verv well with an Refregier fc Ellis 200C| ; [Kaiser, Wilson fc Luppino 2000D ) 



accuracy ~ 10%. We combined difi^erent ray-tracing simula- 
tions with largely difi'erent resolutions (0.1' and 2') to com- 
pare our predictions. The transition from the weakly non- 



should already be able to measure the skewness of the con- 
vergence. 



7 ACKNOWLEDGMENTS 

FB and LvW thank IAS for hospitality, where part of this 
work was done. We thank Bhuvnesh Jain for the invalu- 
able use of his ray-tracing simulations and for useful com- 
ments on the manuscript. RS is supported by endowment 
funds from the Institute for Advanced Study and NSF grant 
PHY-0070928 at IAS. This research was supported in part 
by the Direction de la Recherche du Ministere Frangais de la 
Recherche. The computational means (CRAY-98 and NEC- 
SX5) to do the A'^-body simulations were made available 
to us thanks to the scientific council of the Institut du 
Developpement et des Ressources en Informatique Scien- 
tifique (IDRIS). Numerical computation in this work was 
partly carried out at the TERAPIX data center and on 
MAGIQUE (SGI-02K) at lAP. 



Villumsen, J.V., 1996, MNRAS, 281, 369 
White. W., & Hu, W. 2000, ApJ, 537, 1 

Wittman, D.M., Tyson, J. A., Kirkman, D., DcirAntonio, I., Bern- 
stein, G., 2000, Nature, 405, 143 



REFERENCES 



Bacon, D., Refregier, A., Ellis, R., 2000. 3.stro-ph/0003008 
Bernardeau, P., Valageas, P., 2000, |astro-ph/000627C 



Van Waerbeke, L., Mellier, Y., 1997, A& A, 322, 



Brainerd, T.G., Villumsen, J. V, 



Bernardeau, F 

1 
Blandford, R.D., Saust, A.B 

1991, MNRAS, 251, 600 
Bouchet, F.R., Adam, J.C, Pellat, R., 1985, A&A, 144, 413 
Bouchet, F.R., Juszkicwicz, R., Colombi, S., Pellat, R., 1992, ApJ 

394, L5 
Catelan. P. . K amionkowski 
I ph/OOO5470| 

2000, 



M., Blandford, R., 2000, astro- 



istro-ph/0004151 



Cooray, A., Hu, W., 

Croft, R.A.C., Metzler, C, 2000, fa.stro-ph/0005384 

Erben, T., V an Waerbeke. L.. Bertin. E., Mellier, Y., Schneider, 

P., 2000, [a.stro-ph/0007021 
Fry, J.N., 1984, ApJ, 279, 499 
Goroff, M.H., Grinstein, B., Rey, S.J., Wise, M., 1986, ApJ, 311, 



Hamana, T., Colombi, S., & Suto, Y. 2000, A&A submitted 
Hamilton, A.J.S., Matthews, A., Kumar, P., Lu, E., 1991, ApJ, 
274, LI 



Heavens, A., Refregier, A., Heymans, C, 2000, istro-ph/0005269 

Hui, L., 1999, ApJ, 519, 9 

Jain, B., Seljak, U., 1997, ApJ, 484, 560 

Jain, B., Seljak, U., White, S.D.M., 2000, ApJ, 530, 547 

Kaiser, N., 1992, ApJ, 388, 272 



Kaiser, N., Wilson, G., Luppino, G.A., 2000, istro-ph/0003338 
Miralda-Escude, J., 1991, ApJ, 380, 1 



Moutarde, F., Alimi, J.-M., Bouchet, F.R., Pellat, R., 

A. 1991, ApJ, 382, 377 
Munshi, D., Coles, P., 2000, MNRAS, 313. 148 



: Ramani, 



Munshi, D., Jain, B., 2000a, istro-Dh/9911502 



Munshi, D., Jain, B., 2000b, astro-ph/991233C 



Peacock, J.A., Dodds, S.J., 1996, MNRAS, 280, L19 

Peebles, P.J.E., 1980, The Large-Scale Structure of the Universe 
(Princeton: Princeton University Press) 

Schneider, P., Van Waerbeke, L., Jain, B., Kruse, G., 1998, MN- 
RAS, 296, 873 

Scocciinarro, R., Frieman, J., 1999, A pJ. 520. 35 



astro-ph/0009427 



Scoccimarro, R., Couchman, H., 2000, 

Valageas, P., 2000, A& A, 356, 771 

Van Waerbeke, L., Bernardeau, F., Mellier, Y., 1999, A& A, 342, 

15 
Van Waerbeke, L., Mellier, Y., Erben, T., Cuillandre, J.C, 

Bernardeau, F., Maoli, R., Bertin, E., Mc Cracken, H., Le 

Fevre, O., Fort, B., Dantel-Fort, M., Jain, B., Schneider, P., 

2000, A& A, 358, 30 



